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Abstract 

Entanglement, or quaiitum inseparability, is a crucial resource in quantum information applica- 
tions, and therefore the experimental generation of separated yet entangled systems is of paramount 
importance. Experimeiital demonstrations of inseparability with light are not uncommon, but 
such demonstrations in physically well-separated massive systems, such as distinct gases of atoms, 
are new and present significant challenges and opportunities. Rigorous theoretical criteria are 
needed for demonstrating that given data are sufficient to confirm entanglement. Such criteria 
for experimental data have been derived for the case of continuous-variable systems obeying the 
Heisenberg-Weyl (position- momentum) commutator. To address the question of experimental 
verification more generally, we develop a sufficiency criterion for arbitrary states of two arbitrary 
systems. When applied to the recent study by Julsgaard, Kozhekin, and Polzik [Nature 413, 400 
- 403 (2001)] of spin-state entanglement of two separate, macroscopic samples of atoms, our new 
criterion confirms the presence of spin entanglement. 
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Entanglement, or quantum inseparability, is a profound property of nature that en- 
ables information to be stored, communicated, and processed in a decidedly non-classical 
fashion.[0] Entanglement has long been observed in the states of small numbers of micro- 
scopic objects such as electrons or photons. Only recently have there been efforts to create 
and observe entanglement in the state of massive macroscopic objects, such as the collective 
spins of two separate atomic vapors.|^, § It is important therefore to develop a sufficient 
criterion, which, if satisfied, would unambiguously verify that an experiment has displayed 
entanglement. 

Previous significant work has been done to find a sufficiency criterion that is valid 
for continuous-variable systems obeying the Heisenberg-Weyl (HW) commutator, valid for 
position- momentum variables and, similarly, for light-field amplitudes. 0, ^ Such a con- 
dition is not strictly valid, however, for collective spin systems, although an approximate 
correspondence was proposed for certain special spin states and used to analyse a recent 
experiment by Julsgaard, Kozhekin, and Polzik (JKP).0] This study was aimed at demon- 
strating spin-state entanglement for two separate, macroscopic samples of atoms containing 
around 10^^ atoms each. We derive a sufficiency condition for the existence of entangle- 
ment between two arbitrary quantum systems, including spin systems, in pure or mixed 
states. This allows us, for example, to confirm rigorously the presence of entanglement in 
the experiment of JKP. This new criterion is general, and so may find application in other 
experimental studies. 

Two distinct quantum systems 1 and 2 are said to be entangled if their joint density 
operator p is inseparable, that is, if p cannot be represented as a convex sum of density 
operators pu and p2i for the two physically separated systems, P, 

P = ^PiPli ® P2i, (1) 

i 

with pi a set of non-negative, normalized probabilities. If their joint state is separable (not 
entangled), then it must be possible to express the density operator in the form Eq. (p. 
One physical interpretation of entanglement is that it represents a correlation between two 
systems that is stronger than can exist in any classical (local, realistic) theory.]^ 

A convincing demonstration of entanglement would prove a violation of the separability 
condition Eq. (|l]). In attempting to demonstrate inseparability between the spin variables 
of two separated atomic samples, JKP employ non-local Bell measurements on the spin 
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variables and relate these spin variables to canonical position and momentum operators 
obeying the Heisenberg-Weyl (HW) commutator [qj,Pk] = i^jk {j, k = 1,2). By establishing 
this approximate correspondence, JKP then adapt a criterion by Duan et al. 0] and by 
Simon ||^, which applies to coupled oscillators (and specifically to squeezed light). The 
"HW" criterion that is sufficient for inseparability is 0, |^ 

var(gi + ^2) + var(pi - P2) > 2, (2) 

where var(...) represents the statistical variance. JKP's criterion is an expression analo- 
gous to Eq.(|^), predicated on the assumption that for certain states spin operators can be 
approximately replaced by canonical position and momentum operators. 

Although the shortcut proposed by JKP offers an appeahng connection between criteria 
for demonstrating entanglement in squeezed-hght systems and in spin ensembles, the vahdity 
of this correspondence is far from obvious, and can lead to misconceptions regarding trans- 
formations between different bases that are quite distinct from the Fourier transform nature 
of the canonical position- momentum transformations. Before returning to a consideration 
of entanglement in coUective spin systems, we first establish a criterion for inseparability 
that is applicable to any algebra, including that for spin. We do this by generalizing the 
calculations of Duan et al and of Berry and Sanders |p . 

We consider two systems 1 and 2, and two observables for each, Ai, Bi for system 1 and 
A2,B2 for system 2, that obey [y4i,5j] = 6ijCj . Define linear combinations, 

ii = aAi + (3A2 

v = aBi- (3B2, (3) 

for a, (3 arbitrary real coefficients. Equation ([I|) implies for the variance 

var(w) = ^j9,[a2^(Aii)2), + /52((Ai2)').] 

+S, (4) 

where AAk = — (^fc)p and (...)p denotes an average over p. The quantity S' is S' = 
J^iPiiu)! — (JliPi{u)if' , where (...)» denotes the average over the product density operator 
f>ii ® f>2i- The Schwarz inequality implies in general that S' > 0. Doing the same for v and 
adding the results gives 

var(«) + var({)) > ^^(^.[«^((AiO'). + /^'((Ais)').] + 

pAa\{AB,f). + P\{AB2f)^]), (5) 
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or 

(3'[{{AA,r), + {{AB,r),]. (6) 

Equation is always satisfied for any separable state, with respect to any variables 
(discrete or continuous) belonging to any algebra. If one can measure all the corresponding 
quantities and find a violation of Eq. (^), then one demonstrates that the state is inseparable. 

The general commutator = SijCj imphes the uncertainty relation AAiABi > 

(l/2)Cj, where Cj = \{Ci)\ = \TT{p[Ai, Bi])\, {i = 1,2). This imphes the less restrictive 
relation AAf + ABf > C i, with equality only for AA^ = C» /2. Inserting this into Eq. (||) 
gives, for any separable state, 

var(M) + var({)) > a^Ci + /J^Cs. (7) 

This is our main result. A related criterion has been recently found for the case of pure 
states of spin systems. |Q In the special case a = P = 1, Eq. (|^ gives 

var(ii + ia) + var(5i - B2) > Ci + C2. (8) 

Equation (|^) is not a tight bound. That is, it is necessary for any separable state to 
satisfy Eq. (|^), but it need not be violated for every entangled (i. e., inseparable) state. So 
Eq. (|^ is a necessary but not sufficient condition for separability. A sufficient and necessary 
criterion that is experimentally accessible for spin ensembles is not known. For the special 
case of Gaussian states of Heisenberg-Weyl systems with Ci = C2 = 1, Eq.(||) reduces to Eq. 
(§), which has been shown by Duan et al. and Simon to be a sufficient and necessary 
condition for separability in this case. 

In the JKP study, the variables of interest are the projections J^, Jy, Jz of the coUective 
spins of two atomic samples, 1 and 2. The experiment can be analyzed by choosing 
ii = J,i,fii = J,i,i2 = Jy2.B2 = -J.2. Then C, + C2 = + \{-X2)\ = 2K4i)|, and 

separability requires, from Eq. (P), 

var(J3,i + Jy2) + var(J,i + ^2) > 2|(J^i)|. (9) 

Equation yields a rigorous criterion: if this inequality is violated, then entanglement has 
been demonstrated. 
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This result is similar in form to JKP's Eq. (1), reviewed below, but is distinct in several 
important respects. The first is that our criterion for demonstrating inseparability is ex- 
pressed entirely in terms of the spin operators and does not entail any approximations. This 
result is vahd even without the restriction that a large number of atoms is required. There 
is no recourse, nor any need for recourse, to canonical position and momentum operators 
or to the criterion for squeezed oscillators. The second difference is that Eq. is a vahd 
criterion for arbitrary states not only for certain extremum states as in the criterion of 
JKP. The final difference is that the right-hand side of Eq. @ is the expectation value with 
respect to the state under investigation, rather than being determined by a quantity defined 
in terms of some "classical" value. 

The present result puts on a firm theoretical ground the criterion used by JKP as a 
necessary criterion for separability. The violation of Eq. by the data in the JKP study 
can be taken as an indication of the breakdown of separability. Nevertheless, the question of 
an experimentally accessible, sufficient condition, even for special classes of states (e. g., the 
Gaussian ones for the case of HW systems), is still an open one for the case of spin systems. 

Here we present arguments that one cannot take the approximate correspondence between 
spin variables and HW variables too literally, as it can lead to errors if care is not taken. 
(Our approach avoids the problematic extrapolation of HW results.) For example, large 
errors occur when calculating a change of basis if one uses eigenstates of Jy and as basis 
states and assumes that these transform approximately as HW variables do. These errors 
persist even for the extremum states considered by JKP. 

To review JKP's analysis, the collective-spin vector operator J (total angular momentum) 
of a coUection of N atoms (where may be known or statistically distributed) is defined to 
have x-component Jx = Y^f=i Jx \ similarly for Jy and J z- These components obey the 
algebra [Jy, J z] = iJx, et. cycL, and commute with J^; the number of atoms determines 
the corresponding irreducible representation. For = 2j, one choice for an orthonormal 
basis comprises \j,mx)x which satisfy the eigenvalue relations J'^\j,mx)x = j{j + ^)\j,^x)x 
and Jx\j,fnx)x = fnx\j,fnx)x- For Jx equal to some "large classical" real number (> 0), JKP 
define operators Q = Jy/ \fj^ and P = Jz/ \fjx satisfying [Q, P] = iJx/ Jx |p^ - 

Consider extremum states having narrow support over approximately equal values 
of mx = Jx, where Jx is a large, state-independent real number. Such extremum states 
can be visualized as tightly concentrated near the Jx 'pole' in a space with axes Jx, Jy, Jz, 
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FIG. 1: "Extremum" angular-momentum states having large total J and total Jx=J can be visu- 
alized as occupying the shaded region tightly concentrated near the Jx 'pole' in a space with axes 
Jx, Jy, J z- The quasi-continuous variables Q and P can be thought of as forming approximately 
the Cartesian coordinates of the tangent plane touching the sphere with radius J. 

as illustrated in Fig. |I[ JKP suggest that for such states one can approximate J^/ Jx by 
the unity operator to obtain [Q,P] = i. This commutator, along with Eq.(0), would lead 
directly to the necessary criterion for separabiUty in the form of JKP 's Eq. (1). This result 
is correct in a restricted sense, as noted above. 

Nevertheless, there are difiiculties with taking this approximate approach too Uterally. 
This is evidenced by the fact that a basis transformation between the eigenstates Q = 
Jy/ \fjx. and P = Jz/ \fjx is not given by a Fourier transform, despite the commutator 
between the operators being forced to be a constant, which seemingly imphes that the 
eigenstates have overlap {P\Q) oc exp(— -iPC^). The inappUcabihty of the Fourier transform 
is apparent by attempting this transformation. In the Jy basis, 

oo j 
2j=0 my=-j 
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oo J 

= Y. Y. \j^'^y)yCy{j,my), (10) 

2j=0 my=-j 

and the summation notation means sum j over nonnegative half-integers. For the extremum 
states, with large mean-j value j (say 10^^), the coefficients Cy{j, niy) are non-negligible only 
in the vicinity of = 0. 

The same state represented in the Jz basis is 

oo j 

1^) = E E C7,(j,m,)|j,m,),. (11) 

The two sets of coefficients are related by 

i 

Czii.mz) = E z{j,mz\j,my)yCy{j,my), (12) 

my=-j 

Angular momentum algebra gives z{j,fnz\j,fny)y = '^m^.myl^/^); where the elements of the 



rotation matrix (reduced Wigner function) are 

f^mm'(7r/2) = 



(j + m')!(j -m')! / j+m Wj-mV ^^„^.^^ ^^^^ 



^ (j + m)! (j — m)! ^ + m' — /cy \^ 

The basis transformation Eq. (p!2| ) is entirely different from a Fourier transformation, in which 
the my,mz values would be replaced by quasi-continuous variables niy Jy = Q\fjx^ 
f^z Jz = P\fj'z-, and the transformation would be 

/oo 
dQ{ll2Tif'^eM-iPQ)Cy{Q). (14) 
-oo 

The asymptotic form of the reduced Wigner function Eq.(^) is given in the Appendix, 
where it is seen not to be approximated by the Fourier transform kerneL Furthermore, these 
transformations differ in a quahtative way: Whereas the transformation kernel of Eq. ([1^) is 
necessarily complex, there exists a choice of phase that makes the correct kernel real, as in 
Eq. (H). 

A concrete example, given in the Appendix, illustrates the large errors that can result 
from using the Fourier transform. There we consider a specific state satisfying the as- 
sumed extremum properties (Jy, J z -C Jx)^ which would presumably make the commutator 
[Q,P\ = i approximately correct. Upon making a basis change from the Jy basis to the Jz 
basis, we find, using the correct Eq. (p!2D , that the mean value of Jz is given by a formula 
consistent with J z ^ Jx- However, when (provisionally) using the Fourier transform for the 



basis change calculation we compute a mean value J z = —{t^/2)Jx- This is incorrect, as for 
this state J z must be much smaller than J^.. This demonstrates the complete breakdown of 
a simple, direct replacement of the spin-operator algebra by the HW algebra, leading to the 
need for the more careful derivation we provided in the first half of this paper. 

In conclusion, Eq. (|^) provides a necessary condition for separabihty for arbitrary states 
of two general systems. This condition is accessible to experimental tests in that it involves 
measurements of only several low-order moments. When apphed to collective angular- 
momentum variables in macroscopic atomic systems, the new criterion confirms the one 
used by JKP in their experimental study. [H The problem of finding sufficient conditions for 
special classes of angular-momentum states remains to be solved. 

The complete replacement, for all purposes, of the collective angular-momentum algebra 
by the simpler HW (position-momentum) algebra is not valid, even for extremum states 
that are nearly confined to a small region in angular- momentum space, corresponding to 
highly polarized atomic samples. We do not intend to imply that the use of the approximate 
commutator [Q, P] = i will always lead to large errors. If one evaluates operator moments 
involving only states confined to the proper extremum region, then only small errors are 
incurred, as is well known. We caution, however, that one cannot assume the validity of 
state expansions in basis states having the same properties as Q and P eigenstates. 

Finally, it is interesting to address the question - what states, if any, are conjugate to the 
\j,m)y states through a Fourier transformation? The answer is the SU(2) phase states. In 
the SU(2) phase formulation ||12| one constructs the (2j + l)-dimensional basis from phase 
states, defined as | j, ek)y = (2j + 1)-^/^ J2L=-j e^""^* |j, m)y, with Ok = A;vr/(2j + 1). From this 
we obtain the desired Fourier transform kernel y{j,6\j,'m)y = e^^^/v^^J^FT. Even though 
the basis change from \ j, m) y states to phase states is a (discrete) Fourier transformation, the 
phase operators 0^ constructed for this representation do not naturally yield a commutator 
[Jy, (f)y] = i. Therefore there is not an exact way to use this correspondence to construct an 
equivalent HW algebra. 
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I. APPENDIX 



When j is large and m, m' ^ j the reduced Wigner function Eq. (|I3D is well approximated 
by using Stirling's formula to give 

0(^/2) = ^exp(+|m2 - m"\/2j) cos (^(j + m - m')0 . (15) 

This does not approximate to the Fourier transform kernel. 

As an illustration of the large errors that can arise when using the Fourier transform to 
execute a basis change between Jy and Jz bases, consider the state with 

/ —N j+my j-iriy 

Cy{j,my) = _====_, (16) 
\/ [J + ^yV-U - ^yV- 

where ak = exp(i0fc). The mean values for this state are j = (|aip + |a2p)/2, my = 
Jy = (|aip — |a2p)/2, Jz = \o:ia2\ cos(02 — 0i), and J^ = |aia2| sin(02 — 4>i)- We consider 
states such that 02 — 0i = 7r/2 and — |a2p ^ + |«2p, and hence are in the 
considered extremum class, with Jy, Jz -C Jx- 

In the Jz basis this same state is represented exactly by (using Eqs. (|12|,|13|)) 



c,,,™.) ^ 7'-^'^|;-^-;, (17, 

with Pi = («2 + ai) I \f2 and /^2 = («2 — «i) / \f2. The mean values J^^, J y, J z are unchanged 
by the change of basis, but we now have = J ^ = \aia2\ cos(02 ~ 0i)- 

How does this exact result compare with that obtained by assuming that the HW com- 
mutator is valid, which requires that we transform Eq. (|16]) by the Fourier relation? To 
carry this out we first find an accurate approximation to Eq. (0), using Stirling's formula, 
which gives, for j = j large (e. g. 10^^) and my,my <^ j (and arbitrary phases 0i, 02)? 

Cy{j,my) = C{j){njy^/^exp[-{my - niyf /2]\ exp[i(0i - 02)"^^], (18) 



with C(j) = yexp(— 2j)(2j)2^/(2j)! exp[z(0i + 02)j]5 which is a relatively narrow function 
of j. 

Using Eq. ( p!^ to transform Eq. (plSf ) we find (with the provisional result indicated by 
the tilde) 

exp -i (01 - 02)mj, ,(19) 



1 ^^'^ 



Cz{j,mz) = C{j)[-^ exp 



(mz - 



which means that the mean value of J z is given provisionally by J ^ = riiz = {(f)i — (1)2) J x- 
This predicted value for J z is quite incorrect. As an example consider 02 — 0i = 7r/2. 
Equation (|T7|) predicts (correctly and exactly) that J z has a mean value J z = O, while Eq. 
(p!9|) predicts a mean value J z = Triz = — (7r/2) J^.. This is incorrect, as for this state J z must 
b e much smaller than J^- 

For completeness, the correct Eq. (|17D for the state in the Jz basis can be well approxi- 
mated using Stirling's formula, (since — |/32p -C + |/32p), giving 

Cz{j, mz) = C{j) exp[-(m, - mzf/2j] exp[t{(P[ - <\>'^)m-}^, (20) 

where 0'^ = arg[/?A:]- Equations ( ^O]) and ([T9|) differ in two important ways - both in the 
phase structure and in the predicted mean value of m z- 
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